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Abstract: A promising method for the detection of UHE neutrinos is the Lunar Cherenkov technique, 
which utilises Earth-based radio telescopes to detect the coherent Cherenkov radiation emitted when a 
UHE neutrino interacts in the outer layers of the Moon. The LUNASKA project aims to overcome the 
technological limitations of past experiments to utilise the next generation of radio telescopes in the 
search for these elusive particles. To take advantage of broad-bandwidth data from potentially thousands 
of antennas requires advances in signal processing technology. Here we describe recent developments 
in this field and their application in the search for UHE neutrinos, from a preliminary experiment using 
the first stage of an upgrade to the Australia Telescope Compact Array, to possibilities for fully utilising 
the completed Square Kilometre Array. We also explore a new real time technique for characterising 
ionospheric pulse dispersion which specifically measures ionospheric electron content that is line of sight 
to the moon. 



Introduction 

The origin of the most energetic particles observed 
in nature, the ultra high energy (UHE) cosmic rays 
(CR), which have energies extending up to at least 
2 x 10 20 eV, is currently unknown. Finding the 
origin of these particles will have important astro- 
physical implications. However, direct detection 
of UHE neutrinos is very difficult due to their ex- 
tremely small interaction cross-sections. Instead, 
they may be detected indirectly via observation 
of the Askaryan effect [1] in the lunar regolith. 
Askaryan first predicted coherent Cherenkov emis- 
sion in dielectric media at radio and microwave 
frequencies. Using the Moon as a large volume 
neutrino detector, coherent radio Cherenkov emis- 
sion from neutrino-induced cascades in the lunar 
regolith can be observed with ground based tele- 
scopes. This method was first proposed by Dagke- 
samanskii and Zheleznykh [3] and first applied by 
Hankins, Ekers and O' Sullivan [5] using the Parkes 
radio telescope. 



Coherent Cherenkov radiation is a linearly po- 
larised broadband emission. The spectrum of co- 
herent Cherenkov emission rises approximately 
linearly with frequency until a peak value is 
reached. The peak frequency is determined by de- 
coherence and/or attenuation in the regolith, and 
can vary between a few hundred MHz and approx- 
imately 5 GHz. The dependence of the peak fre- 
quency on shower geometry makes the choice of 
an optimum observation frequency non-trivial [6]. 

Detection Issues 

Lunar Cherenkov emission produces an extremely 
narrow pulse (sub-nanosecond duration). These 
pulses travel through the ionosphere and experi- 
ence a frequency dependent time delay resulting in 
pulse dispersion similar to the dispersion experi- 
enced by pulsar pulses traveling through the inter- 
stellar medium. Therefore lunar Cherenkov pulses 
received on an Earth based radio telescope will be 
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Figure 1 : Pulse dispersed over 50ns (correspond- 
ing to solar maximum) and bandlimited between 
1.2-1.8 GHz. 

dispersed by the ionosphere and further broadened 
by receiver band limiting. 

The ionospheric pulse dispersion must be known in 
real time to maximise the received signal to noise 
ratio and subsequent chances of pulse detection. 
Coherent pulse de-dispersion requires an accurate 
knowledge of the ionospheric dispersion character- 
istic which can be parameterised by the instanta- 
neous ionospheric Total Electron Content (TEC). 
The TEC produces a frequency dependent delay in 
the received signal via Equation 1 

At = 1.34 x 1(T 7 x TEC x z/~ 2 (1) 

where At is in seconds, TEC in electrons per m 2 
and frequency in Hz. 

As Cherenkov radio pulses are much shorter in 
time than any signals normally encountered in ra- 
dio astronomy, real time de-dispersion and detec- 
tion requires broadband technology and innova- 
tions to the current instrumentation. Our experi- 
ment utilises the Australia Telescope Compact Ar- 
ray radio telescope which currently has 600 MHz 
bandwidth at 1.5 GHz but is being upgraded to 
a 2 GHz bandwidth over 1-3 GHz. Due to data 
storage limitations, the only way to exploit these 
new bandwidths is to implement real time de- 
tection algorithms. This requires innovations in 
hardware design to perform signal processing with 
nanosecond timing accuracy. New radio instru- 



ments planned, such as the Square Kilometre Ar- 
ray (SKA), could provide huge advances in col- 
lecting area and technology for the method of lunar 
Cherenkov UHE neutrino detection. 

An array of small dishes is the optimum radio in- 
strument for this experiment. If the dish size is kept 
small enough (~20m at 1.4 GHz), all of the Moon 
will be seen by each antenna's primary beam. In- 
creased sensitivity from using larger dishes is off- 
set by decreased coverage of the Moon due to the 
smaller beam size, and thus a reduced effective 
aperture. Using multiple smaller dishes allows the 
same sensitivity without loss of coverage, provided 
their separation is such that the thermal lunar emis- 
sion — which tends to dominate the system tem- 
perature — is incoherent between antenna. Sig- 
nals can also be added coherently to form multi- 
ple beams around the limb of the Moon, where the 
event rate is maximised, and the array geometry 
can be exploited for RFI discrimination based on 
the signal direction of arrival. 

Current Experiment with the ATCA 

Our current experiment uses the Australia Tele- 
scope Compact Array with a 600 MHz band- 
width. The 600 MHz signal is taken from a main- 
tenance point in the antenna, therefore obtaining 
this bandwidth requires the development of cus- 
tomised detection hardware which must operate 
with nanosecond timing accuracy. Ionospheric 
pulse dispersion must also be corrected in real time 
to maximise the received signal to noise ratio for 
threshold detection [5]. For this stage of our exper- 
iment, we have developed analog de-dispersion fil- 
ters and a field-programmable gate array (FPGA)- 
based trigger system to detect events in real time 
and transmit event data to the site control room via 
an Ethernet link. 

FPGAs are reprogrammable semiconductor de- 
vices, composed of Configurable Logic Blocks 
(CLBs), which can perform simple logic tasks as 
well as more complicated mathematical and signal 
processing algorithms. Their reconfigurable logic 
blocks also provide the ability to optimise hard- 
ware interconnections for each specific algorithm 
implemented. Modern FPGAs can operate at clock 
frequencies of around 500 MHz and have signal 
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processing capabilities that far exceed anything 
previously available for optimised real time opera- 
tion. This high speed operation and optimised pro- 
cessing power make FPGAs an ideal technology to 
perform real time signal processing and pulse de- 
tection. 

The preliminary ATCA experiment makes use of 
an FPGA based sampler board for event trigger- 
ing. Future experiments with the ATCA will make 
use of larger FPGA boards to perform real time 
de-dispersion however for this stage of the experi- 
ment, signal de-dispersion is performed in analog 
microwave filters with a fixed dispersion charac- 
teristic. These filters were designed using a new 
method of planar microwave filter design based 
on inverse scattering [8]. This results in filters 
with a continuously changing profile, in this case 
a microstrip line with continuously varying width. 
The width modulations on the microstrip line pro- 
duce cascades of reflections which sum to pro- 
duce the desired frequency response. For the lu- 
nar Cherenkov detection experiment the desired 
frequency response has a group delay which is 
quadratically chirped according to Equation 1 . As 
the microwave de-dispersion filters have a fixed 
de-dispersion characteristic, an estimate had to be 
obtained for the TEC which would minimise the 
errors introduced by temporal ionospheric fluctua- 
tions. 

Fluctuations in the ionosphere experience a strong 
diurnal cycle and are also dependent on the season 
of the year, phase of the current (1 1-year) solar cy- 
cle and the geometric latitude of observation. Ob- 
servations were planned during night-time hours to 
minimize these fluctuations and as we are currently 
at solar minimum, it was assumed that the electron 
density for our observations could be estimated 
based on measurements from the corresponding 
season last year. This estimate was produced by 
analysis of Vertical TEC (VTEC) data maps which 
were derived from GPS dual frequency signals and 
are available online from NASA's Crustal Dynam- 
ics Data Information System (CDDIS) [7]. 

Our observations were during the nights of May 
5, 6 and 7 2007. These dates were chosen to en- 
sure that the Moon was at high elevation (partic- 
ularly during the night time hours of ionospheric 
stability) and positioned such that we would be 
sensitive to UHE particles from the galactic cen- 



ter. The corresponding CDDIS measurements for 
the month of May 2006 gave an average VTEC of 
7.06 TEC Units (1 TECU - 10 16 electrons per m 2 ) 
over night-time hours (10pm-8am) with a standard 
deviation of 1.3 TECU. The filter design assumed 
a differential delay of 5 ns across the 1.2-1.8 GHz 
bandwidth, which was based on the average VTEC 
value, corrected for slant angle through the iono- 
sphere, and converted to a differential time delay. 
GPS data available post experiment revealed that 
the average VTEC for the nights of our observa- 
tions was actually 7.01 TECU which gave an aver- 
age differential delay of 4.39, with standard devia- 
tion 1 .52, corrected for slant angle. 

Our event trigger has been implemented in an 
FPGA based sampling board which was designed 
as part of the broadband upgrade planned for the 
compact array. The sampler board has two inputs, 
sampled at 2.048 G samples/s, which we are using 
to import two linear orthogonal polarisations of the 
600 MHz RF signal. The sampled inputs are then 
multiplexed into parallel 512 MHz data streams for 
input into a Xilinx FPGA which performs pulse de- 
tection. Detection is performed via thresholding 
on both polarisation streams. The sampled data is 
buffered for 2 microseconds so that data surround- 
ing candidate events can be sent back and recorded 
in the case of a detection. 

Event data is transferred to a central processing site 
via a Gbit Ethernet connection for off-line process- 
ing. The de-dipsersion filters and CABB sampler 
boards were installed on three antennas so that co- 
incidence testing and direction of arrival discrimi- 
nation could be performed during off-line process- 
ing. As each antenna triggers independently, trig- 
ger timing information and nanosecond synchro- 
nisation is essential for both of these processing 
stages. 

Our current sampler boards do not have the capa- 
bility of time stamping with nanosecond accuracy, 
and so alternative methods of time stamping had 
to be investigated. As the system temperature is 
dominated by thermal emission from the Moon, 
the array operates in the intensity interferometer 
regime [2], and there is a small amount of corre- 
lation between the signal received at each antenna. 
This level of correlation can be exploited to de- 
termine relative timing between event buffers sent 
back from different antennas. 



Developments in Nanosecond Pulse Detection 



Future Experiments and Developments 
in Pulse De-Dispersion 

Future improvements to the ATCA UHE neutrino 
detection experiment include using 5 antenna, per- 
forming coherent signal combination in real time 
to enable coincidence testing and an increase to 2 
GHz bandwidth. Real time coherent signal addi- 
tion will be possible as part of an upgrade planned 
for the ATCA which will include the installation of 
powerful new FPGA based back-end receiver hard- 
ware. This hardware can also be used to implement 
real time de-dispersion algorithms and we have de- 
veloped a technique for obtaining measurements 
of the ionospheric TEC which are both instanta- 
neous and line-of-sight to the lunar observations. 
The ionospheric TEC can be deduced from Fara- 
day Rotation measurements of a polarised source 
combined with geomagnetic field models, which 
are more stable than ionospheric models [4]. The 
Faraday Rotation induced in a radio wave is related 
to the ionospheric electron content via Equation 2 

ft = 2.36 x 10 4 z/- 2 / N(s)B(s) cos Ods (2) 

J path 

where il is the rotation angle in radians, v is the 
signal frequency in Hz, TV is the electron density 
m~ 3 , B is the geomagnetic field strength in T, 
is the angle between the direction of propagation 
and the magnetic field vector and ds is a path el- 
ement in m. Our novel approach is to use this 
technique with the polarised thermal radio emis- 
sion from the lunar limb as our polarised source, to 
obtain instantaneous, line-of-sight TEC measure- 
ments. This makes the lunar Cherenkov technique 
extremely attractive for UHE cosmic ray and neu- 
trino astronomy as it removes the need for search- 
ing in dispersion space. 

The next phase in our experiment involves the 
Square Kilometre Array (SKA). The SKA is a new 
generation radio telescope array which will be one 
hundred times more sensitive than the best present 
day instruments. Current designs proposed for 
the SKA consist of large numbers ( 10 4 ) of small 
dishes (6- 12m) to achieve a square km of collect- 
ing area in the 0.1-3 GHz range which is critical 
for UHE neutrino experiments. To gain the advan- 
tage of the large number of small dishes offered by 



the proposed SKA designs, signals have to be com- 
bined and analysed with nanosecond timing accu- 
racy. This will involve forming phased array beams 
in real time using special purpose beam forming 
hardware. We will require enough beams to cover 
the visible surface of the Moon over the entire fre- 
quency range, and will thus have increasing resolu- 
tion with frequency. Using these methods, we ex- 
pect the SKA sensitivity to reach the level at which 
a flux of UHE neutrinos will be detected. 
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